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[57] ABSTRACT 

A non-invasive in-vivo method of analyzing a bone for 
fracture risk includes obtaining data from the bone such 
as by computed tomography or projection imaging 
which data represents a measure of bone material char- 
acteristics such as bone mineral density. The distribu- 
tion of the bone material characteristics is used to gener- 
ate a Finite element method (FEM) mesh from which 
load capability of the bone can be determined. In deter- 
mining load capability, the bone is mathematically 
“compressed”, and stress, strain force, force/area ver- 
sus bone material characteristics are determined. 

9 Claims, 4 Drawing Sheets 
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METHOD FOR IMPROVED PREDICTION OF 
BONE FRACTURE RISK USING BONE MINERAL 
DENSITY IN STRUCTURAL ANALYSIS 

5 

This invention was made with Government support 
under Grant No. 1 P01-DK-39964 awarded by the Na- 
tional Institute of Health and Subcontract 956084 
awarded by NASA-JPL. The Government has certain 
rights in this invention 1 ° 

BACKGROUND OF THE INVENTION 

This invention relates generally to the analysis of 
bone structure for determining risk of fracture, and 
more particularly the invention relates to an in vivo 15 
non-invasive method of determining bone fracture anal- 
ysis by obtaining a measure of bone mineral density, 
establishing bone structure, and analyzing the structure 
for load carrying capability. ^ 

DESCRIPTION OF PRIOR ART 

The amount of bone, or bone mineral density (BMD), 
or other properties of bone such as the speed of sound, 
are used conventionally to evaluate the skeletal status of 
individuals in an effort to predict the susceptibility of 2 
that bone, or by inference other bones in the patient, to 
fracture under minimal applied loads such as those en- 
countered in normal daily living. This evaluation is 
done in the common use by passing a collimated beam ^ 
of radiation through an object such as a person or inani- 
mate object and measuring the transmitted photons on 
the opposite side of the object. The intensity of trans- 
mitted photons is compared to the intensity transmitted 
through a known object for the purpose of calibration. 35 
The calibrated intensity is used to describe some prop- 
erty of the object, such as bone mineral in the path of 
the beam. Alternatively, scattered or reflected radiation 
can be analyzed instead of transmitted radiation (such as 
Compton scattered photons or ultrasound). Many dif- 40 
ferent types of apparatus can be used to do this, includ- 
ing but not limited to: computed tomography scanners, 
x-ray or radioisotope source projection imaging systems 
(including film), single-beam scanners or ultrasound 
devices. For all devices, the eventual outcome is a quan- 45 
titative measurement of an average property of the bone 
measured. 

Some devices generate an image of the object proper- 
ties measure. Such images can be analyzed using various 
methods to give a regional distribution of the object 50 
properties, and these distributions are at times com- 
pared to distributions derived from population- 
averaged data or to data derived from prior studies of 
the same object. Such comparisons are done in a region- 
by-region or point-by-point basis but are not used to 55 
derive specifically the distribution of material proper- 
ties. The values derived from such analyses are com- 
pared to a large, previously-derived database of values 
from other individual measurements, both from individ- 
uals with fractured bones (osteoporotics) and those 60 
without fracture (normals) to determine how well the 
value measured can classify a given patient as osteopo- 
rotic or normal. In some cases, combinations of mea- 
surements from different regions or different bones are 
used to try to improve this classification procedure. 65 
However, in virtually all cases the use of such measure- 
ments explains only 50-60% of the variance in the pre- 
dictive capability for classification. 
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In order to try to improve this predictive capability, 
researchers have attempted to correct the measured 
bone mineral density by some parameter of geometry, 
such as cross-sectional area. Other investigators have 
developed analytic models of bones in an attempt to 
incorporate structure information into the analysis of 
bone strength. While these models can provide some 
additional information, they require that each bone be 
modeled directly, and thus they cannot be used in a 
practical setting. The best data regarding strength of 
bones comes from direct mechanical testing of bone. It 
is impossible to test bones directly in the living subject, 
with the exception of long bones such as the ulna or 
tibia where bending stiffness properties have been deter- 
mined. No other researchers have proposed a general 
method for direct determination of the mechanical 
strength or failure properties of individual bones in 
individual subjects. 

SUMMARY OF THE INVENTION 

An object of this invention is to provide a method to 
measure non-invasively (in vivo) the strength of an 
individual bone in an individual patient. This measure- 
ment can then be used to determine whether or not the 
bone will fracture under specified loading conditions 
such as those normally seen in daily living. It can also be 
used to estimate fracture risks under abnormal loading 
conditions such as occur in falling, jumping or during 
athletic events or heavy training regimens. 

The invention uses the distribution of physical prop- 
erties of bone measured non-invasively in an individual 
and mathematical analysis of that distribution to predict 
the risk that a bone may fracture under applied loads. 
The use of such methods relates to the clinical disease of 
osteoporosis, or in general metabolic bone diseases, 
although by inference such methods can also be used to 
evaluate bones in any situation where the amount of 
bone may be compromised, such as bone metastases in 
cancer, multiple myeloma, or Paget’s disease. In a pri- 
mary application, 3-dimensional quantitative computed 
tomography data acquired using a conventional CT 
scanner are used to determine the distribution of bone 
mineral density, this distribution is used to define bone 
material properties, and the finite element method of 
analysis is used to determine structural properties of the 
whole or a part of the bone. Other applications include 
the use of any 2-D or 3-D noninvasive method to deter- 
mine the distribution of bone material properties, such 
distribution measured in an individual then analyzed by 
the finite element method to predict structural proper- 
ties. 

In the general application, the invention relates to the 
use of non-invasive methods to determine material 
properties of an object followed by use of these proper- 
ties as input to a finite element method (FEM) of analy- 
sis to determine failure modes. In some cases, the rela- 
tionship between the measured property (e.g. density) 
and the actual property of the object input into the 
analysis (e.g. Young’s Modulus) may require informa- 
tion acquired under other, conditions. 

The invention and objects and features thereof will be 
more readily apparent from the following detailed de- 
scription and appended claims when taken with the 
drawings. 
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BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 is a schematic illustrating one apparatus to 
acquire data for the distribution of material properties in 
accordance with one embodiment of the invention. 

FIG. 2 illustrates one method of practicing the inven- 
tion. 

FIG. 3 A and 3B illustrate one means of obtaining 
bone fracture determination using the invention. 

FIGS. 4A-4C illustrate the mathematical analysis of 10 
a vertebra in accordance with the invention. 

FIG. 5 is a schematic illustrating another apparatus 
for acquiring data for the distribution of material prop- 
erties in accordance with another embodiment of the 
invention. 

DETAILED DESCRIPTION OF ILLUSTRATIVE 
EMBODIMENTS 

An advantage of the invention is its ability to use 
specific information obtained about the distribution of 20 
materia] properties in a bone in a patient to determine 
noninvasively the strength of that bone. This informa- 
tion is then used to predict risk of fracture under speci- 
fied loading conditions. Specifically, the distribution of 
bone material properties determined non-invasively 25 
using one of a number of techniques (computed tomog- 
raphy, projection imaging, ultrasound, magnetic reso- 
nance) is used as input to a finite element analysis of 
structural strength, and other parameters such as load- 
ing conditions and boundary conditions are also in- 30 
eluded in the model as needed. Using mathematical 
methods contained in commercially-available or spe- 
cially written computer programs, the model of a bone 
can be incrementally loaded until failure, and the yield 
strength determined. Alternatively, other mechanical 35 
properties can be measured this way. This invention is 
the first method to incorporate the distribution of bone 
material properties of any individual bone into an analy- 
sis of that object’s mechanical characteristics. In order 
to practice this invention, it is required to measure the 40 
distribution of a parameter of bone, relate this parame- 
ter to material properties, generate a matrix containing 
geometrical and material properties of the object, and 
subject this matrix to defined loading conditions using 
the known finite element method of analysis to deter- 45 
mine mechanical characteristics of the object as a 
whole. 

A computed tomography (CT) imaging system as 
employed in accordance with one embodiment of the 
invention is illustrated in FIG. 1. A radiation source is 50 
mounted on a conventional gantry with a radiation 
detector on the opposing side of the patient. The CT 
imaging system is used to obtain an image of the patient 
within the reconstruction area, such image representing 
an x-y map of the x-ray attenuation properties of the 55 
patient and having a finite thickness in the z-axis defined 
as the “slice thickness”. The CT scanner table is incre- 
mented by a preset amount, normally the slice thickness 
and a second image is obtained. The process is repeated 
until enough images have been obtained so that a whole 60 
bone or a desired part of a bone is included in the data 
set. 

The output of the apparatus of FIG. 1 is convention- 
ally a single number or several numbers representing an 
average material property of the bone. Our invention 65 
takes the data from such apparatus prior to the reduc- 
tion to average properties, converts the data to a matrix 
of the geometrical distribution of material properties, 
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uses the matrix as input to conventional finite element 
analysis software, and uses the finite element method 
(FEM) of analysis to generate mechanical properties. 
See for example Huiskes “A Survey of Finite Element 
5 Analysis in Orthopedic Biomechanics: The First De- 
cade,” Pergamon Press Ltd., 1983 pages 385 - 409 . Such 
properties are output from the finite element analysis 
directly or by analysis of preliminary output such as a 
stress vs. strain curve. A property such as the yield 
stress of the bone is then used to characterize the pa- 
tient’s risk of fracturing the bone, or by inference, of 
fracturing other, similar, bones. The CT scanner is cali- 
brated so that the x-ray attenuation properties can be 
related to properties of the bone, such as bone mineral 
15 density. Individual regions of bone mineral density can 
be measured from this data set (Cann CE, Genant HK: 
“Precise measurement of vertebral mineral content 
using computed tomography,” J. Comput. Assist. To- 
mogr. 4 : 493 - 500 , 1980 ). 

In FIG. 2, the 3-D CT data set is acquired and the 
information about the distribution of bone is preserved. 
The original CT data are then processed using known 
methods to convert each point in the distribution to 
bone material characteristics (BMC), and the bone can 
be separated from surrounding materials not used in the 
analysis, if desired. The bone data set is rotated in 3 
dimensions to an orientation required for the mathemat- 
ical loading analysis. The coordinates of each distribu- 
tion CT point and their material properties are used to 
generate a conductivity matrix or mesh describing the 
exact relationships among the regions of the bone. The 
mesh is input to finite element method (FEM) analysis 
software using appropriate boundary conditions. FEM 
models have heretofore been used in orthopedic biome- 
chanics to model bones for use in designing orthopedic 
prostheses. The bone is mathematically “compressed” 
for structure analysis. Force versus displacement (stress 
vs. strain) information about the elements in the mesh 
(that is, the regions of the bone) is determined and such 
information can be plotted on a graph, for example in 
FIG. 3, or analyzed on the computer to determine the 
yield strength of the bone (the point at which elastic 
deformation converts to plastic deformation). In FIG. 
3A stress vs. strain in normal bone is plotted, and in 
FIG. 3B stress vs. strain in osteoporotic bone is plotted. 

This is further illustrated in FIGS. 4A, 4B, and 4C. 
FIG. 4A illustrates the original scan plane orientation 
for CT scans of vertebrae. The bone material character- 
istics of a vertebra of interest are determined and ro- 
tated in orientation for uniform loading as shown in 
FIG. 4B. The elements of the rotated vertebra after the 
mathematical compression test are illustrated in FIG. 
4C. 

FIG. 5 illustrates analysis of a hip bone (proximal 
femur) using a two dimensional projection raster scan in 
accordance with another embodiment of the invention. 
Again, the scan data for the bone is extracted and con- 
verted to bone material characteristics (e.g. bone min- 
eral density). Other known techniques can be used for 
obtaining the initial data representing the bone charac- 
teristics including DP A, DEXA, and QDR. 

There has been described an improved method of 
prediction of bone fracture risk by obtaining data from 
which distributed bone material characteristics such as 
bone mineral density can be determined and from which 
structural analysis can be made. The method is espe- 
cially useful in diagnosis and treatment of patients with 
or predisposed to osteoporosis, but the method can be 
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used for other individuals such as athletes and astro- 
nauts in training. Further, the FEM model can be modi- 
fied to reflect aging of an individual. 

While the invention has been described with refer- 
ence to specific embodiments, the description is illustra- 5 
tive of the invention and is not to be construed as limit- 
ing the invention. Various modifications and applica- 
tions may occur to those skilled in the art without de- 
parting from the true spirit and scope of the invention as 
defined by the appended claims. 

What is claimed is: 

1. A method of determining fracture risk of a bone in 
vivo under normal loading conditions based on the 
density and geometric distribution of bone mineral com- 
prising the steps of 

non-invasively obtaining x-ray attenuation data from 
said bone representing a measure of bone mineral, 

converting said data into a bone mineral distribution 
characteristic, 

generating a finite element method model based on 
said bone mineral distribution characteristic, and 

determining load capability of said bone using said 
FEM model. 

2. The method as defined by claim 1 wherein said step 
of non-invasively obtaining x-ray attenuation data in- 
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eludes acquiring computed tomography (CT) data for 
said bone. 

3. The method as defined by claim 1 wherein said step 
of non-invasively obtaining x-ray attenuation data in- 
cludes obtaining projection image data. 

4. The method as defined by claim 2 wherein said step 
of obtaining x-ray attenuation projection image data 
includes a raster scan of a particle beam through said 
bone. 

5. The method as defined by claim 1 wherein said step 
of determining load capability of said bone includes 
mathematically compressing said FEM model. 

6. The method as defined by claim 5 wherein said step 
of determining load capability further includes deter- 

15 mining stress and strain versus bone material character- 
istics for said bone. 

7. The method as defined by claim 6 wherein said step 
of determining load capability further includes deter- 
mining failure load for said bone. 

20 8. The method as defined by claim 7 and further 

including the step of altering said FEM model to repre- 
sent aging of said bone. 

9. The method as defined by claim 1 and further 
including the step of altering said FEM model to repre- 

25 sent aging of said bone. 
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